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ABSTRACT 

The intense irradiation received by hot Jupiters suppresses convection in the outer layers of their 
atmospheres and lowers their cooling rates. "Inflated" hot Jupiters, i.e., those with anomalously large 
transit radii, require additional sources of heat or suppressed cooling. We consider the effect of forced 
turbulent mixing in the radiative layer, which could be driven by atmospheric circulation or by another 
mechanism. Due to stable stratification in the atmosphere, forced turbulence drives a downward flux 
of heat. Weak turbulent mixing slows the cooling rate by this process, as if the planet was irradiated 
more intensely. Stronger turbulent mixing buries heat into the convective interior, provided the 
turbulence extends to the radiative-convective boundary. This inflates the planet until a balance 
is reached between the heat buried into and radiated from the interior. We also include the direct 
injection of heat due to the dissipation of turbulence or other effects. Such heating is already known 
to slow planetary cooling. We find that dissipation also enhances heat burial from mixing by lowering 
the threshold for turbulent mixing to drive heat into the interior. Strong turbulent mixing of heavy 
molecular species such as TiO may be necessary to explain stratospheric thermal inversions. We show 
that the amount of mixing required to loft TiO may overinflate the planet by our mechanism. This 
possible refutation of the TiO hypothesis deserves further study. Our inflation mechanism requires 
a deep stratified layer that only exists when the absorbed stellar flux greatly exceeds the intrinsic 
emitted flux. Thus it would be less effective for more luminous brown dwarfs and for longer period 
gas giants, including Jupiter and Saturn. 

Subject headings: diffusion - opacity - planet-star interactions - planets and satellites: atmospheres 
- radiative transfer - turbulence 



1. INTRODUCTION 

Hot Jupiters — giant planets receiving intense irradi- 
ation from their host stars — are the best characterized 
class of exoplanets. Their proximity to their host star 
yields frequent transits of and occultations by their host 
star that are observed over a range of wavelengths. Two 
remarkable features of hot Jupiters are their inflated radii 
and the variety of infrared emission signatures, some of 
which have been interpreted to reveal stratospheric ther- 
mal inversions. 

Many hot Jupiters have larger radii than standard cool- 
ing models predict, even with the intense irradiation from 
the host stars included in the radiative transfer. This 
implies a mechanism that injects heat into and/or re- 
tards the loss of hea t from the planets' interiors. See 
iFortnev et all jjOOg) for a review of proposed mecha- 
nisms; 

iGuillot k Showman! (|2002l hereafter GS02) argued 
that atmospheric winds, driven by intense irradiation, 
could explain inflated radii. In their model, the kinetic 
energy of the winds dissipates as heat below the pen- 
etration depth of starlight. However the energy need 
not be deposited into the convective interior (a common 
misconception). Dissipating energy in outer radiative 
layers suffices to delay planetary contraction. Turbu- 
lence, which the winds can trigger via Kelvin-Hclmholz 



instabil ities, is an efficient me chanism to dissipate kinetic 
energy (jLi k GoodmaiJlMtOh . MHD drag is an alterna- 
tive dissipation mechanism provided weather-layer winds 
extend to the high-pressure, metallic zone of hydrogen 
(jPerna et alJl2010l : iBatvgin k Stevensorjl201fl ). 

Thermal inversions, i.e. regions where the atmospheric 
temperature rises with height, may also implicate tur- 
bulent mixing in hot Jupiter atmospheres. Transit 
spectra of several hot Jupiter s have been interpreted 
as being thermally inverted (iRichardson et al.l 120071 : 
iBurrows et"aLl 120071 : iKnutson et al.ll20I0D Th ese obser- 
vation s appear to confirm the predicion of iHubenv et al.l 
(2003) that molecular absorbers, mainly TiO, in the 
stratosphere could generate inversion layers. However 
vapor phase TiO could rain out of the upper atmosphere 
if it condenses in cold traps. Turbulent mixing can coun- 
t eract this settl i ng. 

iSpieeel et all (|2009l hereafter S09) showed that eddy 
diffusion coefficient^] of K zz « 10 7 — 10 11 cm 2 /s are 
needed to maintain sufficient stratospheric TiO for ther- 
mal inversions. The range of K zz values accounts for 
the varying extent of cold traps in planets with different 
thermal profiles and the size of grains that condense. S09 

1 Eddy diffusion models small-scale turbulent processes in anal- 
ogy to molecular diffusion, with tracers fluxed down their mean 
gradients. 
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argued that the need for strong mixing renders the TiO 
hypothesis "problematic," pending improved estimates 
of K zz . One goal of our study is to determine if large 
K zz values are energetically problematic. 

Sulfur has also been proposed as a high-altitude ab- 
sorber. Photochem ical models of sulfur abundances 
(|Zahnle et al.l 120091 ) also include turbulent mixing at a 
strength K zz w 10 7 cm 2 /s, though the dependance on 
K zz is unclear. Eddy diffusion is also used in brown 
dwarf models to explain disequilibrium chemical abun - 
dances ([Griffith fc Yelldll99l iHubenv fc Burrowsll2007l ). 

Turbulence not only mixes chemical species, it also 
transports heat. This paper develops a model that in- 
cludes this turbulent heat flux in the radiative layers of 
hot Jupiters. While convection drives an outward flux 
of energy, forced turbulence in stably stratified regions 
drives a downward flux of energy. This effect is distinct 
from — though it accompanies — the dissipation of tur- 
bulence as heat, which we also include. 

By altering the flow of energy, we change the cooling 
and contraction rates of hot Jupiters. For modest lev- 
els of turbulent diffusion, the outward radiative flux is 
partially offset by the downward flux of mechanical en- 
ergy. This reduces the net cooling flux from the convec- 
tive interior, which self-consistently pushes the radiative 
convective boundary (RCB) to higher pressure^ 

For sufficiently strong eddy diffusion, the downward 
flux of energy exceeds the outward radiative flux that 
a planet of fixed entropy can provide. In this case the 
turbulent heat flux flows into the convective interior, in- 
creasing the internal entropy and inflating the planet. A 
schematic of this mechanism is shown in Fig. [TJ Because 
higher entropy planets are more intrinsically luminous, 
inflation leads to an equilibrium between turbulent heat 
burial and radiative losses. 

Our mechanism bears some resemblance to the run- 
away greenhouse. In the latter, the atmosphere is com- 
posed of a greenhouse gas in vapor pressure equilibrium 
with a large, surface volatile reservoir. The cooling emis- 
sion to space emanates from a pressure ~ g/n, with sur- 
face gravity g and (Rosseland mean) opacity k. The 
emission is independent of the surface temperature for 
optically thick atmospheres. Vapor pressure equilibrium 
determines the temperature at the emission level, thus 
limiting the cooling radiation that the atm osphere can 
achieve (Kombavash^ Il967t UngersolU 119690 . If the ab- 
sorbed sunlight exceeds this limiting cooling emission, 
the surface temperature increases until either the volatile 
reservoir is depleted or the atmosphere becomes suffi- 
ciently transparent to the surface blackbody emission. 
The role of limiting cooling flux in the runaway green- 
house is played in our mechanism by the cooling flux 
of the core. The role of absorbed sunlight is played by 
the downward, mechanical flux of energy. If the latter 
exceeds the former, the planet heats up by increasing 
the core entropy until energy balance can be achieved. 
The analogy is somewhat incomplete, however, in that 
the traditional runaway greenhouse involves a radiative- 
thermodynamic feedback which does not exist in our 
mechanism. 

2 For simplicity we will describe convectively stable regions as 
"radiative," even when we include the transport of heat by both 
turbulence and radiation. 
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Figure 1. Schematic of the mechanical greenhouse effect to in- 
flate hot Jupiters. A downward flux of heat (large black arrow) is 
driven by turbulence in the convectively stable "mixing layer" and 
deposited in the deep interior. This downward flux can balance 
or even exceed the convective losses (gray overturning arrows). 
Atmospheric circulation ("winds") launched near the photosphere 
drive turbulence in the mixing layer. Other mechanisms, such as 
non-linear gravity wave interactions, could also drive the turbulent 
flux. 

This paper is organized as follows. In §Yl we review 
standard radiative equilibrium models. We add turbu- 
lent heat transport and energy injection to our model in 
§31 We derive our mixing length formulae for the tur- 
bulent transport of heat in ij3.il and present the model 
equations and our solution methods in §3.21 We provide 
a prescription relating turbulent diffusion and dissipa- 
tion in stratified atmospheres in §3.31 We present and 
analyze our model results in §31 We first treat constant 
diffusion (partly to connect with S09) and ignore energy 
dissipation in §33] We then add complexity by consider- 
ing a spatially varying K zz in £14.21 and including energy 
dissipation in M4.31 We discuss consequences of chang- 
ing the opacity law in i|4.4l We compare our results to 
dynamical simulations of atmospheric circulation in i j5.ll 
and to the TiO diffusion needed for thermal inversions 
in £|5.2[ We summarize our results and their implications 
in §H 

2. STANDARD ATMOSPHERIC MODELS 

We start with a review of the standard radiative 
transfer approximations used in this work (§23]) and 
apply them to radiative equilibrium solutions f i|2.2[) . 
We will introduce our n otation and parameter choices. 
lArras fc Bildstenl ([20061 hereafter AB06) present a simi- 
lar analytic model, which they compare to global models 
with detailed opacities and equation of state (EOS). 

2.1. Radiative Transfer 

Our goal is to understand energy balance. We focus 
on the deep atmosphere which is optically thick both to 
incoming stellar irradiation and the planet's emitted flux. 
Here the equation of radiative diffusion 

dT _ F lad 
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relates the outgoing radiative flux F la d to the variation of 
temperature T with pressure P via the radiative diffusion 
coefficient 

P) 



^rad 



where cr is the Stcphan-Boltzmann constant. Hydrostatic 
balance, dP/dz — —pg allows pressure to replace height z 
as the vertical corrdinate, with p the atmospheric density. 
We hold gravity g constant, invoking the plane-parallel 
approximation for thin atmospheres. 

For the Rosseland mean opacity, our calculations will 
use a power law, 



K = K 1 P a T fi = K — — 



P 



T V 



T 2k J 



(3) 



The two forms are equivalent, with the constant m being 
more compact, while k has units of opacity and is nor- 
malized to Pkb = 1 kbar and T%k — 2000 K. Unless stated 
otherwise, calculations will use n oc P, i.e. a — = 
as a rough approximation to collision induced molecular 
opacity. Our normalization choice of k q = 0.18 cm 2 /g 
will be justi fied below f H2.2p . We discuss alternate opac- 
ity laws in £14.41 While realistic opacities are only well 
approximated by power laws over a limited range, this 
considerable simplification is useful for developing intu- 
ition. 

At the top of our atmosphere, we set the temperature 
to Tdeem an ap proach used in AB06 and advocated by 
llro et all (|2005f ). This approach is valid when the inci- 
dent stellar flux exceeds the emitted radiation, resulting 
in a deep isothermal region at the top of the optically 
thick atmosphere. In this physical situation, the pre- 
cise location of the upper boundary is not important, 
as we explain furt her in ||2.2| W e point the reader to 
IHansenl (|2008l ) and lGuillotI (|2010D for sophisticated ana- 
lytic treatments of radiative equilibrium, which are very 
useful for interpreting co mputational models. Note that 
some (jBaraffe et al.ll2003h but not all (jSeager fc Sasselovl 
I2000D detailed radiative transfer solutions show an ex- 
tended isotherm in hot Jupiter atmospheres. One differ- 
ence between models is the choice of non-gray opacities, 
which affects the depths at which starlight in different 
frequency intervals is absorbed. 

The incident stellar flux averaged over the full plane- 
tary surface, Fi TI = <jT^ gives a characteristic tempera- 
ture 



TI w 2000 K 



-1/4 



M l/6pl/3 
- iV2 *,©- r day 



(4) 



where stellar mass and luminosity, M* © and T*,© are 
normalized to solar values, and the orbital period Pday 
is normalized to a (24 hour) day. Horizontal tem- 
perature gradients are important for driving winds in 
the weather layer. However these winds efficiently 
smooth temperature gradients at pressures > bar, where 
timesc ales for advection are shorter than for radiative 
losses (Showm an et al.l I2009D . Thus ID models are ap- 
propriate for basic considerations of energy balance. 

Because the planet is not a perfect blackbody, Td CC p 
may not match . Greenhouse or anti-greenhouse effects 
depend on the relative transparency of the atmosphere 
to stellar and emitted longwave radiation. (To be clear, 



we are now referring to standard radiative effects, not 
the mechanical greenhouse.) If incoming starlight pen- 
etrates below the infrared photosphere, then the green- 
house effect gives Td eop > T*. If, however, significant 
incoming radiation is absorbed above the photosphere, 
a stratospheric ther mal inversion gives Td ee p < T*. Sec 
iHubenv et alj (|2003[ ) for a more quantitative analysis. In 
most of our examples we adopt Tdccp = 1500 K, as appro- 
priate for short (~ I day) orbital periods with a thermal 
inversion, or for longer periods with no thermal inver- 
sion, a greenhouse effect and/or a more luminous host 
star. 

Giant planets, including hot Jupiters, become unstable 
to convection at depth. The lapse rate of the atmosphere 



dlnT 3kPF t 



ad 



d\nP 16p crP 4 



(5) 



characterizes its stability, and the final equality follows 
from equation (fT) . Convection occurs where V > V a d- 
We set V a d = 2/7, the adiabatic index of an ideal di- 
atomic gas. In reality, non-ideal interactions lower V a d 
at the high pressures of exoplanet atmospheres, and pro- 
mote convection. At the top the atmosphere, where the 
optical depth t = nP/g k, 1 and F ra d <C crT^ eep , equa- 
tion (|5|) shows that V <C 1 and the atmosphere is in- 
deed stable and nearly isothermal. For reasonable opac- 
ity choices, V increases with depth, and gives a transition 
to convection (even under the ideal gas approximation). 

In convective regions we set V = V a d, i.e. an adiabatic 
profile with T oc P Vad . The efficiency of convective en- 
ergy transport makes the modest super-adiabaticity neg- 
ligible. The level of the adiabat is determined by the in- 
ternal entropy. A global calculation of entropy is beyond 
our ill ustrative scope. Instead, motivated by IHubbardI 
(1977), we label our adiabats by Ti, the temperature it 
would have at Pi = 1 bar pressure, even though the adia- 
bat likely does not extend to such low pressure. We define 
a reference entropy (per unit mass) S^ef, corresponding 
to Ti = 250 K. Relative entropy values for different T\ 
are then computed as 



AS = S - SVcf = C P ln(Ti/250 K) 



(6) 



where the specific heat (at constant pressure) Cp = 
i?/V a d is assumed constant. For the gas constant R — 
/cs/(/xm p ) we use a mean molecular weight p = 2.34 
times the proton mass m p . 

The stable atmosphere matches smoothly onto the 
convective adiabat at the radiative-convective boundary 
(hereafter RCB). Since the temperature T c and pressure 
P c at the RCB lie on the interior adaibat we require 



' 'ft 



(7) 



The location of the RCB is crucial for global evolu- 
tion. The secular cooling of the convective interior is 
determined by the radiative flux, F c , leaving the RCB. 
Combining equations d3j, ([5]) and then (|Tj) at the RCB 
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gives 



Fc. 



Va d 16ffaT c 4 - 



3kiP< 

V ad 16gcr 

3«i 



1+a 



Ti 



4-/3 



P 4 



pVadpV 



P 



6/7 



(8) 



(9) 



with 



(l + a)/(4-/3) = l/2. 



(10) 

Core flux increases with the interior entropyO Pushing 
the RCB to higher pressures decreases the core flux if 
Voo > V a d- This condition is satisfied for our opac- 
ity choice, and is generally required for a transition to 
convection (as we show shortly). We emphasize that the 
dependance of P c on P c is independent of the mechanism 
that changes P c , though previous works have mostly con- 
sidered irradiation. These basic considerations are use- 
ful in interpreting numerical studi es of planetary cool- 
ing histories and r adii evolution (iBurrows et al.l 120001 : 
IBaraffe et al.M2003t IChabrier et al.ll200l ~ 

2.2. Radiative Equilibrium Solutions 

We now apply the two approximations of radiative 
equilibrium (RE) to the stable layer. First, radiation 
is the only relevant energy transport mechanism. Thus 
P ra d in equation (P) is the total flux of energy. Sec- 
ond the flux is constant through the radiative layer with 
Pad = Fc, the flux from the convective interior. This 
assumes that local (thermal) energy release is negligible. 

Fig- HI plots radiative equilbrium atmospheres for k oc 
P, with two values of Pdccp matched on to interior adia- 
bats labeled by T\. We obtain analytic RE solutions by 
integrating equation (JTJ) with P ra< j = F c and T — Pdccp 
at P = to get 



T = T< 



deep 



' ad 



Voo ~ V 



ad 



1+a 



1/(4-/3) 



(11) 



This solution uses equation ([8j and we have imposed the 
requirement P(P C ) = T c to find 



Fc Pleep 



v~ 



1/(4-/3) 



(12) 



, Voo — V a d , 

A valid solution — one that transitions to convection — 
thus requires Voo > V a d and a > — 1 (which together 
assure j3 < 4). As Fig. [2] shows, F c increases with Pd CC p 
but is independendent of the interior entropy. 

The RCB sinks to larger pressure as entropy decreases 
or as Tdeep increases, 

l/Vad 



P c = fc v Pi 



T, 



deep 



(13) 



which follows from equations ([7]) and (TT2"j) with the con- 
stant 



Voo/[V ad (l+a)] 



2.1 



s Voo — V a d , 

3 The numerical scaling in equation Q ignores the effect that 
higher entropy would lower gravity by inflating the planet. This 
effect cancels when computing the total luminosity, which is ulti- 
mately more important. 
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Figure 2. Radiative equilibrium (RE) atmospheres with deep 
isotherms of T^ eBp = 1500 K(blue curves) and 2000 K (dotted green 
curves) matched onto internal adiabats (dashed red curves) with 
entropy increasing from bottom to top. Grey dots mark the loca- 
tion (T c and P c ) of the radiative-convective boundary (RCB), and 
squares show Pdeepi where V = V a( j/2. 



The core flux for RE atmospheres follows from equa- 
tions (151). OH) and (ITB1 as 



P = k F - 



Ti 



(l+a)/Vad 



nl-V ad /V 



deep 



deep 



(14) 



This gives the well known result that increased irradia- 
tion reduces the cooling of the planet, while higher en- 
topy planets are more luminous. The constant 



hp* = 



16ctV 



ac 1 



3P 



1+a 



1 



Vad 
Voo 



Voo/V a 



Equation (| 14[) is consistent with, but more specific than, 
equation © in assuming that RE sets the location of P c . 

We chose the parameters for Fig. [2] — used throughout 
this work — by roughly matching the analytic solutions 
to more detailed hot Jupiter models, as in AB06. 

We constrain the entropy parameter F\ by appealing 
to the typical P c ~ 1 kbar location of the RCB in hot 
Jupitcrs with modest, i.e. Jovian, entropies. With T\ 



260 K, we reproduce a 1 kbar RCB for F^e 



1500 



K. We also consider larger values of T\ to describe more 
inflated planets, but keep P c ^> 1 bar. 

The normalization of the opacity determines the core 
fluxQ Requiring P c = er(100 K) 4 for the standard pa- 
rameters and P c = 1 kbar, gives k q = 0.18 cm 2 /g for 
g = 10 3 cm 2 /s. We emphasize that this is not a realistic 
opacity law (in particular it is too low at small pressures) . 
We are merely choosing parameters that allow the simple 
analytic model to mimic properties of more detailed hot 
Jupiter models. 

The lapse rate for RE solutions is (from eq. [Tl]) 



V = V 



ad 



(P/Pc) 1 - 



1 



V, 



Va 



(P/Pc) 



1+a 



(15) 



4 Remarkably, k does not affect the location of the RCB along 
a given adiabat, only T^j-p and the power laws are required. Over 
long times though the opacity at the RCB affects entropy evolution 
and thereby RCB location. 
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demonstrating that V = V a d at P = P c and that the 
solution becomes isothermal V — > at low pressures. 
Smooth opacity laws give a monotonic increase in V with 
P. Opacity windows give more complicated profiles of V, 
including multiple zones of convection (see M4.4[) . 

We define Pdeep, the effective depth of the isothermal 
layer, as the location where V = V a d/2. This occurs at 



P 



deep 



V 



ad 



2V„o - V 



1/(1+") 



P c « 0.55P C 



(16) 



Our definition of Pdeep differs from AB06, who define 
Pdeep as a characteristic scale that might exceed P c . 
We now revisit the validity of applying the boundary 



condition T 



T, 



at P = 0. Due to the isother- 



mal layer at low pressures, applying the boundary con- 
dition at any P <C Pdeep gives indistinguishable solu- 
tions. However solutions are only physically valid in 
optically thick regions, for P ^> Pthick ~ g/ftmin ~ 
10 mbar[ft m j n / (0.1 cm 2 /g)] _1 . The relevant opacity K m - m 
is the smaller of the opacities to starlight and emit- 
ted radiation near Pthick- Indeed the penetration of 
starlight below the infrared p hotosphere c an push the 
top of the isotherm to ~ 1 bar (|GuiHotil2010D . As long as 
Pthick <C Pieep, solutions are physically consistent below 

Pthick- 

3. ENERGETICS OF TURBULENT RADIATIVE 
LAYERS 

We now generalize the radiative equilibrium model to 
include two effects. First we allow for turbulent eddies 
to drive an advective heat flux P e ddy The total flux 



P — Prad + P 



eddy • 



(17) 



includes the radiative and eddy contributions. Second, 
we allow for the release of energy at a rate e. In steady 
state this heating is balanced by cooling from the diver- 
gence of the total flux, 



dF 
dP 



(18) 



Sources of e include the viscous dissipation of turbulence 
f £|4.3[) . brea king waves (|Showman et a l. 2009) and ohmic 
dissipation (|Batvgin &: Stevensonll2010T ). 

While equation ([T]) still describes the radiative flux, 
P ra d is no longer constant with height. The fractional 
contribution of P ra d to the total flux can vary, and the 
total flux itself can vary. To proceed further we require 
a model for P c ddy and e. 

3.1. Turbulent Heat Transport 

We derive P e ddy using basic elements of mixing length 
theory. This theory is usually applied to convectively un- 
stable regions, but instead we apply it to forced turbu- 
lence in convectively stable regions. We will show that in 
this case the energy flux is inwards. We leave the forcing 
mechanism of turbulent motions unspecified, and their 
strength a free parameter. 

We consider parcels of gas which conserve entropy and 
maintain pressure equilibrium with their surroundings as 
they exchange position over a vertical distance t and then 
dissolve. These parcels contain excess heat Sq = pCpST 



ST= — 



f dT 




dT\ 


\ dz 


ad 


dz J 


IT dS 




c P 


dz 





with: 



(19) 

where Sf gives the difference of any quantity / between 
the parcel and its surroundings. For stable stratifica- 
tion (dS/dz > 0) rising parcels (£ > 0) cool and sinking 
parcels heat. We express the heat flux, P G ddy = wSq with 
w the characteristic eddy speed, in terms of the turbulent 
diffusion, K zz — w£, as 

Pcddy = ~K zz pT- 



dz 



= -K zzP g[l- — 

Vad 



(20) 



The flux is always negative for stable stratification. It 
vanishes at the RCB where dS/dz — (and V = V a d). 
We do not model overshoot, which could allow energy 
exchange (in either direction) between convectively sta- 
ble and unstable zones. In the upper isothermal regions 
Pcddy ~K ZZ P, declines in magnitude with height, un- 
less K zz increases with height to compensate, as we will 
consider in t j4.2l 

The above assumption that parcels conserve entropy 
assumes a radiative cooling time longer than the eddy 
turnover time. We are ignoring radiative losses during 
eddy motion in this work, because we consider optically 
thick regions with long cooling times. This assumption 
will eventually break down in optically thin regions, no- 
tably the inversion layer itself, which is strongly strat- 
ified. Whether eddy fluxes affect the structure of the 
inversion layer is left to future work. 

To understand the energetics of Pcddy, we analyze its 
divergence, which describes cooling and (when negative) 
heating 

n cPS 
dz 2 



dF, 



eddy 

dz 



99 dS 
R dz 



dK zz T dS 
dz dz 



(21) 



The terms on the right hand side represent cooling rates, 
— Sq, which we relate to rates of work, Sw, using the first 
law, Sq = 5e—5w. Thus — Sq = Sib/V^d since the internal 
energy, Se = pCyST = (1 — V a d)<5<7 with Cy = Cp — 
R. The first term in equation (|2ip arises from buoyant 
work Swb = Spgw = K zz pg{dS/dz)/Cp, where Sp/p — 
—8T/T by pressure equilibrium. The first term on the 
RHS of equation (PH]) is thus — SqB = <5?«B/V a d. This 
buoyant cooling will be evident in the stratified regions 
(P< -^clccp ) of Fig. 1 

The second term represents the tendency of mixing 
to heat by filling in entropy minima (for d 2 S/dz 2 > 
0). More specifically it arises from compressional work, 
which vanishes for a constant entropy gradient because 
the work done on rising and sinking parcels cancels. For 
varying dS/dz consider two parcels that arrive at z, one 
from above and one from below. Compressional work is 
done on the parcels at a rate 



PV -v± = - 



P5 P ± 
p St 



Pi dS 



W 
1 



(22) 



Cp dz 

where the top (bottom) sign refers to sinking (rising) 
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parcels, X7 -v is a velocity divergence, and St = £/w gives 
the expansion rate. The net work is the sum of these 
terms, Swq = —V a dK zz pTd 2 S/dz 2 . We identify the sec- 
ond RHS term in equation (|21j) as — Sqc = <5wc/V a d. 
This compressional heating dominates in deeper regions 
(P > Pdeep) of Fig. [51 The third and final term repre- 
sents a flux imbalance that arises from non-uniform eddy 
diffusion as in £|4.2I 



3.2. Model Equations and Self- Similar Solution 
Technique 

Our atmospheric model is described by equations ([TJ), 
([!?)) ■ (fTS)) and (l2Tfl) which reduce to the following pair of 
coupled ODEs 



F + F is 



dT__ 

dP ~ fc rad + F iso P/(V ad T) 
dF _ e 
dP~~~g 



where 



:P9 



(23) 
(24) 

(25) 



is (minus one times) the isothermal limit of the eddy 
flux. The thermal profile in equation (1231) describes the 
combined effects of radiative and eddy fluxes. Ignoring 
mixing (Fi SO — ^ 0) recovers standard radiative diffusion 
of equation (fT]) . Strong mixing (F[ so — > oo) creates an 
isentropic profile (V — > V a( j). 

Solving the coupled ODEs for T and F requires pre- 
scriptions for K zz and e and boundary conditions. As 
with the RE solutions, we fix Td ee p at the top of the at- 
mosphere, and match onto an adiabat with a given T± at 
the bottom. This matching generally requires iterative 
techniques. We avoid this complication by finding self- 
similar solutions. This technique is only possible because 
of the idealizations — notably power law opacities and 
ideal gas EOS — described in iJ2.ll 

To find self-similar solutions, we normalize all quanti- 
ties to their value at the RCB. As with the analytic RE 
solutions, we do not know where the RCB is located un- 
til we find the solution. We integrate outwards from the 
RCB with trivial boundary conditions: the dimensionless 
T and F are unity. 

We set the strength of diffusion not with a physical 
value for K zz , but by the parameter 



4> c 



K zz p c g 
F c 



(26) 



the ratio of Fi so to F = F ra d at the RCB. For dissipa- 
tion we must similarly specify eP/(gF) at the RCB. See 
appendix [Bl for details. 

To get a physical solution, we scale a dimensionless 
solution to any desired value of Td ee p and T\. Setting 
T = Td ccp at P = 0, the solution gives T c at the RCB. 
Specification of T\ then fixes P c via equation ([7]) . We 
then determine F c and K zz via equations © and (p?6"l) . 

3.3. Relating Diffusion and Dissipation 

Turbulence that gives rise to diffusion, K zz , will also 
dissipate at some rate e. We now consider a prescription 
that sets a lower bound on e from turbulence. We also 



allow for stronger heating, perhaps from non-turbulent 
sources. 

In a Kolmogorov cascade, the dissipation rate e = w / £ 
and the diffusion K zz = vol give a simple relation 
e = K zz /t 2 , with t Q = l/w the turnover time of the 
dominant eddies. Unfortunately we lack a reliable model 
for eddy timescales. Moreover turbulence in a stratified 
atmosphere is likely anisotropic and not well described 
by Kolmogorov scalings. Fortunately, our results are 
not very sensitive to this limitation, as we will show in 
subsequent sections. Eddies with long turnover times, 
t Q > 1/A organize into horizontally extended pancakes, 
where the squared buoyancy frequency is 



n2 = rt [Vad ~ v] 



(27) 



Assuming that the buoyancy frequency sets the rele- 
vant timescale, t a = 1/A, gives a dissipation rate 

ebuoy « K ZZ N 2 . (28) 
In strongly stratified, isothermal regions the buoyancy is 



Adccp — g/ \/CpT dccp ~ c s / H 



(29) 



with c s the sound speed and H the scale height. For 
sub-sonic turbulence with w <C c s , our prescription gives 
£ ~ w/N -€i H . This is consistent with the expectation 
that stratification limits the vertical extent of turbulent 
structures to less than a scale height. Forced turbulence 
with t < 1/N is also possible. Since this would give 
even stronger dissipation, we consider ebuoy a reasonable 
lower bound on dissipation for stratified turbulence. 

Near the RCB, as A — > it is unreasonable to expect 
the dissipation to vanish entirely. Thus we also include 
a floor to the dissipation 



feK 



zz Jy deop 



(30) 



with the dimensionless normalization / c giving the ratio 
of e G to ebuoy m isothermal regions. If we use rotation 
as the other relevant timescale, then the floor would be 



quite low with f t 



n 2 /N 2 

/ deep 



I0" 3 P 



-7/3 



prescription considers both terms, 

^ T" ebuoy ■ 



day 



Our full 



(31) 



as discussed in 



4. RESULTS 



4.1. Constant Diffusion, No Dissipation 

We describe solutions to the atmopheric model of £j3j 
We start by considering turbulent mixing with a con- 
stant diffusivity K zz and no dissipation, i.e. e = 0. This 
requires integration of equation ([23]). While the decay of 
turbulence always gives some dissipation, the effect on 
energetics can be small (as we will show in H4.3[) . 

Our main goal is to constrain the amount of turbulent 
diffusion that can be maintained in convectively stable 
regions. The simplest and most generous — i.e. allow- 
ing the largest levels of turbulent diffusion — constraint 
comes when we ignore dissipation. 

An appeal to energetic balance shows why this should 
be the case. Recalling that turbulence drives a downward 
eddy flux, F e ddy < 0, in convectively stable regions, we 
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rearrange flux balance as — -Feddy = F Ta ,d — F. For a 
large downward eddy flux we need the total flux F to be 
small compared to the outgoing radiative flux. Ignoring 
dissipation helps in this regard by preventing F from 
increasing through the layer. Pushing the RCB to high 
pressure also helps by lowering the constant F = F c from 
the RCB. The solutions below show that strong mixing 
does indeed push the RCB to high pressure. 

4.1.1. An Upper Limit to K zz 

Fig. [3] shows the effect of varying K zz while hold- 
ing irradiation (Td oe p = 1500 K) and internal entropy 
(Tj = 250 K) hxed. As K zz increases, the downwelling 
eddy flux pushes the RCB to higher pressures, P c , and 
lowers the flux from the interior, F c . These effects are 

coupled since T c cx P c 6 ^ 7 along a fixed adiabat (eq. [5]). 
Turbulent mixing — like strong stellar irradiation — re- 
duces the planet's cooling rate. 



100 
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20000 
15000 




0.01 



Kzz [cm 2 /s] 

Figure 3. As K zz increases the RCB moves to high pressure (P c , 
dashed green curve) while the core flux (F c , blue curve) drops. 
Both diverge at at finite K zz = K zz<CI it ~ 1660 cm 2 /s, when 
dissipation is ignored (e = 0). This upper limit varies with in- 
ternal entropy and and T decp as shown in Fig. [4] Quantities are 
plotted relative to the K zz = case where P c = 1.1 kbar and 
F c = ct(100 K) 4 . 




Figure 5. Profiles of a stirred atmosphere (blue curves) with 
K zz ss K zz , cr it and no dissipation compared to the RE case with 
K zz = (dotted black curves). Both join an adiabat with T\ = 250 
K (dashed red and gray curves). Mixing pushes the RCB (gray 
dots) to high pressure (which remains finite because K zz is 0.1% 
below K ZZ}CI it). (Left:) The temperature profile shows modest 
changes near Pdeep (blue squares), as shown in the inset. (Right:) 
Lapse rate V, relative to V a( j- Turbulent diffusion smoothes the 
transition towards the adaibat. The inset plot of 1 — V/V a( j shows 
the marginal stability of the stirred atmosphere up to high pres- 



AS [k B /m p ] 
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250 350 500 750 
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Figure 4. Maximum eddy diffusion in the stable layer (i^2 2)C rit)i 
vs. internal entropy for T^ccp = 1000, 1500 and 2000 K (dotted red, 
blue and dashed green curves, respectively). See equation ( 1 3 2 D for 
analytic fits. Entropy is given in terms of Ti [eq. [7] and relative to 
the Ti = 250 K reference (top axis). 



At a critical value of K zz = K zz>a a, P c diverges to 
infinity while F c drops to zero. This upper limit to dif- 
fusion is -?Czz,crit ~ 1660 cm 2 /s for the adiabat and Td CC p 
chosen in Fig. [3l Of course a real planet cannot extend 
to infinite pressure, to say nothing of the plane-parallel 
approximation. The point is that our steady state model 
cannot energetically support K zz > K zz cr i t . 

Stronger turbulence could in principle exist, since we 
are saying nothing about what forces turbulence. In a 
non-equilibrium state with K zz > K zz ^ cr i t , the down- 
welling eddy flux would then increase the internal en- 
tropy, and inflate the planet. 

Fig. 2] shows that higher entropy planets have a higher 
steady state i^ zz , C rit- Thus by inflating the planet, strong 
mixing brings the planet's energy balance into equilib- 
rium. An ultimate upper limit to is that the planet 
not over-inflate and exceed its observed radius. The K zz 
values invoked in S09, from 10 7 to > 10 10 cm 2 /s would 
imply significant or (on the upper end) excessive infla- 
tion. For comparison AB06 showed (see their Fig. 11) 
that entropy changes of AS ~ k^/m v (the scale on our 
Fig. [4]) can expand a hot Jupiter's radius by ~ 10 — 25%. 
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Accurate determination of the maximum K zz allowed for 
a given planet requires more detailed modeling (including 
global structure with realistic opacities and EOS) than 
we perform. However our results strongly suggest that 
K zz values invoked in the literature have significant, or 
even excessive, effects on energetics. 

Fig. 2] also shows that if ZZjCr it increases with decreas- 
ing Tdcop- Thus our constraints on mixing are much more 
stringent for hot Jupiters than for more distant plan- 
ets, including Jupiter itself. Recall that thermal inver- 
sions lower Tdeep and that mixing can sustain thermal 
inversions by keeping opacity sources aloft in the strato- 
sphere. A planet can accommodate strong mixing with 
some combination of thermal inversions to lower Tdeep 
and increased internal entropy. It is hard to predict if 
thermally inverted planets should be more inflated - 
due to the presumed presence of turbulence — or less 
inflated — because lower Td ee p promotes cooling and in- 
hibits our mechanical greenhouse effect. Observations do 
not indicate an obvious correlation. Planets with signa- 
tur es of inversions exh ibit varying degrees of inflation. 
See IMiller et al.l (2009) for a comparison of observed to 
model radii of transiting planets. 

While strong mixing pushes P c — > oo, the depth of 
the isothermal layer, Pdeep, is relatively unchanged by 
mixing. (We explore structure in detail below.) Thus 
planets with higher entropy or lower Td CO p have shallower 

isothermal layers. Specifically Pd CO p oc (Tdeep/Pi) 1 ^*" 1 
from equations (fT5|) and (|TB| . 

It is hardly surprising that planets which can ac- 
commodate more mixing (larger K ZZtCr n) have shallower 
stratified layers to mix (smaller Pdeep)- In principle 
strong mixing could destroy the deep isotherm altogether 
by pushing Pdeep to optically thin regions. This proba- 
bly requires unrealistically large core entropies. Alterna- 
tively as interior temperatures rise, blackbody emission 
at depth may find opacity windows at wavelengths longer 
than the infrared. 

4.1.2. Structure and Energetics of Stirred Atmospheres 
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Figure 6. Energetic balance for the profiles in Fig. [5] In radiative 
equilbrium {black dotted curve) the radiative flux is constant. With 
eddy diffusion and no dissipation the total flux {gray dotted curve) 
is constant. The radiative flux {blue curve) leaving the interior is 
low, but -F ra[ j first increases and then decreases outward, peaking 
near Pdeep {blue square). The downwclling eddy flux {dashed red 
curve) offsets changes to the radiative flux. 



We now consider the structure and energetic balance of 
"stirred atmospheres" with K zz w K zz ,ciit- We compare 
these solutions to standard RE atmospheres of £|2.2l with 
K zz = 0. Since our model is self-similar, the behavior 
is independent of the parameters (Ti, Td CC p) chosen for 
illustration. 

Fig. [5] {top panel) shows that the temperature profile 
of the stirred atmosphere is very similar to the RE case. 
In the stirred atmosphere, the RCB lies at much higher 
pressure below an extended "pseudo-adiabat, "0 which 
lies very close to the original adiabat. The inset to Fig. [5] 
{top panel) focuses on the region near Pdeep (which drops 
slightly to 480 bar from the original 610 bar) where the 
stirred atmosphere is hotter, by at most 60 K. The stirred 
atmosphere is slightly colder below 250 bar, though this 
difference of at most 3 K is not visible. The stirred at- 
mosphere is very modestly thicker, by 0.06Pdecp, where 

Pdeep = PPdccp/,9 ~ 0.008Pj. 

Fig. [5] {bottom panel) plots the lapse rate. Mixing 
smoothes the transition towards the adiabat. The in- 
set shows the smooth decline of 1 — V/V a d along the 
pseudo-adiabat, which gradually reduces the amplitude 
of the downwelling eddy flux, from equation (|20p. 

Fig. [6] shows that the energetics of the stirred at- 
mosphere differs significantly from the RE case. With 
K zz — there is no eddy flux and the radiative flux 
is constant down to the RCB, here at P c w 1.1 kbar. 
The stirred atmosphere has a deeper RCB which reduces 
the core flux significantly. We could push P c — > oo and 
P c — > with a 0.1% increase of K zz all the way to 
Kzz, cr it, but choose not to for visualization. 

The radiative and eddy fluxes change with height. 
Their sum — the total flux — remains constant be- 
cause we ignore dissipation. Fig. [6] shows that the 
fluxes behave differently above and below Pdeep j i- e - 
along the isothermal and pseudo-adiabatic regions, re- 
spectively. Along the pseudo-adiabat, the radiative flux 
declines with depth as P ra d oc P~ 6 / 7 , as we derived 
for the core flux in equation©. The radiative cooling 
{dF Ta a/dP < 0) in this region balances heating by eddy 
diffusion (dP c ddy/dP > 0). Achieving this energetic bal- 
ance requires only modest changes to the T-P profile. 
Since P e ddy scales as 1 — V/V a d, it is very sensitive to 
small changes in V along the pseudo-adiabat (see eq. [2"0] 
and the bottom inset of Fig. [5]) . 

The energy balance along the isotherm, i.e. above Pdeep 
is different. With V <C 1, the eddy flux P dd y ~ 
—pgK zz fx — P grows in magnitude with depth (a dif- 
ferent scaling holds if we vary K zz with height, see 34. 2p . 
This localized cooling (dP c ddy / dP < 0) balances radia- 
tive heating (dP ra d/dP > 0). The decline in radiative 
flux with height again requires only modest changes to 
the thermal profile. From equation ([5]) , P ra d is sensitive 
to small changes in V <C 1. 

We can now simply estimate K ZZjCI it using our knowl- 
edge that P c — >• and that mixing only modestly changes 
the RE profile . The transition region near Pdeep is cru- 
cial. Here the eddy flux reaches its peak negative value 
Peddy,decp ~ - Pdcc P gK zz /2. The thermal profile con- 
strains P la d to be roughly F C {K ZZ = 0), the core flux of 
the RE atmosphere. We set P ra d,dee P ~ 2F C {K ZZ = 0) 

5 This is not to be confused with the pseudo-adiabat that de- 
scribes moist convection in Earth's atmosphere. 
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to account for the slightly hotter atmosphere near Pdeep- 
Energetic balance, P e ddy,dee P + Prad,dccp = P c -> 0, then 
gives 

\ 2+a 



K zz . 



crit 



AF C {K ZZ = 0) 

Pdeepff 



1/V a 



1/V^-l/V'^ 



T 



(32) 



deep 



where VJ^ = (2 + a)/(5 — /3) and our parameters give 
results of Fig. 0J 



K zz , cr it oc if /T^^ . These scalings agree with the 



4.2. Spatially Varying K zz 

The above analysis f t|4.1[) shows that upper limits on 
a constant K zz are set by the balance of eddy and ra- 
diative fluxes near Pdecp, which itself scales with internal 
entropy and Tdeep- To test the robustness of this find- 
ing, we include a depth dependence to K zz oc P*». For 
winds driven near the photosphere, i.e. the top of our 
atmospheres, one might expect stronger diffusion in the 
upper atmosphere, i.e. £ < 0. On the other hand if tur- 
bulence is triggered by shear layers with the convective 
interior, perhaps £ > 0. As discussed in Sj6] [and section 
5?] detailed dynamical simulations can help determine 
plausible diffusion profiles. 

Fig. [7] shows the effect of varying £ with other param- 
eters fixed (at our standard values of Td CC p = 1500 K, 
T 7 ! = 250 K, a = 1, £ = as in e.g. Fig.©. These plots 
show the strongest possible mixing, which (as we found 
for constant K zz ) drives the RCB to infinite depths and 
reduces the core flux to zero. 

The maximum mixing near Pdeep is relatively un- 
changed, except when the mixing at the top of the at- 
mosphere is quite strong. Quantitatively we compare 
values of K zz ,deep, defined as the maximum value of K zz 
at a reference P = 550 bar, which is Pdeep of the radia- 
tive equilibrium atmosphere. For constant K zz we found 
K ZZlClit = K zzAccp = 1665 cm 2 /s. For mixing that in- 
creases with depth as £ = 0.5, 1.0, and 1.5, if zz ,deep de- 
clines by a modest 5%, 6% and 5%, respectively. When 
mixing declines with depth as £ = —0.5, —1.0, and —1.4, 
if 22 , deep increases by 14%, 58% and 300%, respectively. 
We cannot consider models with £ < —1.5 because they 
do not approach an adiabat at depth. The bottom panel 
of Fig. [7] shows that the approach to the adiabat is al- 
ready quite gradual for £ = —1.4. We explore this issue 
further in appendix [A] 

The top panel of Fig. [7] shows the flux profiles for sev- 
eral £ values. The plot shows both radiative and eddy 
fluxes, which obey P ra d = — Peddy because the net flux, 
P — > when mixing pushes the RCB to infinite depth. 
We also plot the radiative flux for the reference radiative 
equilibrium model (horizontal black dotted line) without 
any mixing. The explanation for these flux profiles mir- 
rors the discussion of Fig. [5] in £14.11 At high pressures the 
flux is controlled by the radiative flux along the adiabat. 
Increasing or decreasing the mixing with depth has little 
effect on the deep eddy flux. Changes to K zz are com- 
pensated by (1 — V/V a d) — see equation (|20|) — which 
is small and sensitive to slight changes in V close to the 
adiabat. 

The flux in the low pressure, isothermal region scales 
as — Pcddy oc pK zz cx P 1+ i (see eq. [20]). This explains 
why the flux components, P ra d = —Peddy, increase with 
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Figure 7. Models with a depth-dependent X zz oc P - ^ for £ = 
1,0,-1 and —lA(dotted green, dashed red, dot-dashed purple and 
solid blue curves, repspectively) and no dissipation. Strong mixing 
pushes P c — > oo while Pdecp ls plotted with colored squares. (Top:) 
Radiative flux, which is equal and opposite to the eddy flux. Lower 
f values give strong mixing and larger fluxes at the top of the 
atmosphere. (Bottom:) Thermal profiles show that strong upper 
atmosphere mixing (low f ) heats the upper atmosphere and results 
in a more gradual approach to the adiabat. 



height if £ < —1. Driving larger radiative fluxes in the 
upper atmosphere requires a steeper dT/dP. The tem- 
perature profiles in Fig. [7] (bottom panel) reflect this. 
The £ = — 1 and especially the top £ = —1.4 curves 
are noticeably hotter at intermediate pressures and have 
smaller Pdeep- We thus find that mixing at the top of 
the atmosphere is more effective — compared to uniform 
or bottom-focused mixing — at lifting (i.e. heating) the 
T — P profile. The additional heat in this case is provided 
by a downward flux of mechanical energy across the top 
boundary. 

4.2.1. Limits on Mixing Near the Photosphere 

We now consider what might constrain K zz near the 
top of the atmosphere, since we find that internal entropy 
mostly constrains diffusion near Pdecp- Our £ < — 1 solu- 
tions in Fig.[7]show that strong mixing at the top requires 
a large flux of mechanical energy at the top of the atmo- 
sphere. Weather-layer winds are a plausible source of 
mechanical energy, and they are generated by the atmo- 
spheric heat engine driven by insolation. The thermody- 
namic efficiency of all planetary atmospheres in the Solar 
System is of order 1%. If we restrict the magnitude of 
Peddy to & fraction /* ~ 1% of the insolation P* ~ a Tj ccp 
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scaled for a downwelling flux that originates at a 0.1 bar 
photosphere. While the efficiency /* and mechanisms 
of generating a downward mechanical flux are uncertain, 
the energetic difficulties of mixing at K zz ;§> 10 9 cm 2 /s 
is evident. 

Alternatively, we could attempt to constraint mixing in 
the upper atmosphere by appealing to our £ = — 1.4 solu- 
tion. This gives the largest downward eddy flux (Fig. 
top). Smaller £ values would be needed for a larger 
flux, but these do no give consistent solutions. In ap- 
pendix [A] we analyze the £ ss —1.4 limit and argue that 
it may not be physically significant. Ignoring this con- 
cern, the diffusion near the top of our £ = —1.4 solutions 
is K zz w 5 • 10 7 (P/bar) -14 cm 2 /s. A larger internal en- 
tropy could support a larger K zz (as we showed for con- 
stant K zz models in Fig. HJ. Therefore this constraint 
is not inconsistent with equation ()33[) , which is a more 
robust constraint. 
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Figure 8. Increasing dissipation reduces -R" ZZjC riti shown for fixed 
T dcop = 1500 K and T\ = 250 K. The solid blue curve only includes 
a constant floor to the dissipation, e , while the dashed curve also 
includes our prescription for dissipation in stratified regions. 



4.3. Including Dissipation 

We now consider the effect of adding dissipation to our 
models with eddy diffusion. The total flux F will now in- 
crease with height due to dissipation. The coupled equa- 
tions ([23"f and (pM)) govern the steady state structure. To 
understand the effect of dissipation on the turbulent heat 
flux, we return to the simpler case of spatially uniform 
K zz . 

With dissipation we still find an upper limit to diffu- 
sion, 2£" zz ,orit, f° r a given Ti and T dccp . However K zz , crit 
declines with increasing dissipation. Fig. [5] shows this for 
both a dissipation rate e Q that is constant with height 
(solid blue curve) and the full dissipation prescription 
(dashed red curve, see eq. |31|). The full prescription 
includes ebuoy, our estimate of the minimum dissipation 
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Figure 9. Similar to Fig. [6] except dissipation is included. The 
total flux is no longer constant and the core flux is larger. Both 
of these effects reduce the magnitude of F e[ jdy which can no 
longer offset as much of F ra( j. Consequently -K" ZZlC rit is reduced 
to 900 cm 2 /s. The dissipation profile e = e + e^uoy includes a 
floor e a = 5 X 10 -5 erg/(gs), that corresponds to f e = 0.01, i.e. 
weaker dissipation near the RCB than in stratified regions. 

due to stratified turbulence. This additional dissipation 
further reduces K ZZtCI i t . 

We confirm that the dissipation-free estimates of 
K zz , a it m previous sections represent a conservative up- 
per bound. Lowering if Z2jC rit means that weaker tur- 
bulent diffusion will inflate the planet. Admittedly, the 
cases shown in Fig. [8] do not prove that all dissipation 
profiles will lower K zz clit . However we investigated the 
effects of both spatially- varying K zz (as in tj4.2p and also 
different profiles of e. In all cases adding dissipation re- 
duced K ZZ CI it from the dissipation-free value. 

4.3.1. How Dissipation Lowers K Z z,crit 

We explore the energetics of how dissipation lowers 
Kzz,cmt- Fig. [S] shows flux balance with dissipation, and 
can be compared to Fig. [6l Notice that the peak value 
of — -Fcddy now falls well short of i^ad at the relevant 
pressure, Pdoop- This is because the total flux 
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now includes the integrated dissipation, causing F to 
greatly exceed F c , the small loss of heat from the core. 
From the RCB to Pdcep, -Frad also increases with height. 
This rise is not significantly affected by dissipation, with 

Fmd oc P c 6 ^ 7 along the pseudo-adiabat as before. The 
magnitude of — F e ddy — ^rad — F is smaller at Pdeep 
because dissipation increases F without comparably in- 
creasing F rad . 

We find that Pdccp (and also pdeep) do not significantly 
change when we add dissipation. Indeed the T — P or 
V profiles for the solution in Fig. [5] are indistinguishable 
from the stirred atmospheres in Fig. [SJ except the RCB 
is not pushed as deep, "only" to P c f» 11 kbar. 

The limiting value of K zz = 2|P cd d y |/(Pdcc P g) drops 
because the smaller eddy flux is not compensated by a 
lower pdeop- It seems possible that some dissipation pro- 
file could heat the atmosphere and lower Pdeep and pdeep 
enough to increase if ZZjC rit- We did not find this to be the 
case. One reason is that too much dissipation can affect 
the location of the RCB, a subject we address below. 
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For completeness we explain the flux balance at small 
pressure illustrated in Fig. [9] The decline in — F ddy ~ 
pgK zz with height in isothermal regions again results 
from declining density. However -F ra d does not decline 
towards low pressure (as was seen in Fig. [5]), because the 
total flux F is larger with dissipation. The fact that the 
escaping flux matches the flux from the radiative equi- 
librium solution [dotted black line) is a coincidence (with 
some significance, see below). This coincidence occurs 
because the dominant dissipation is e ~ ebuoy Larger 
or smaller choices of dissipation would give a larger or 
smaller (respectively) net F and escaping F ra d- 

This coincidence has some significance. Downward 
eddy fluxes reduce the loss of heat from the core, as we 
have discussed extensively. However we now see that this 
loss is matched by the flux due to the dissipation of that 
turbulence — provided our prescription for the minimum 
£buoy is correct. This replacement is intriguing, but does 
not alter our discussions of evolutionary consequences: 
turbulent dissipation in radiative regions is powered not 
by interior heat, but by external means (such as forced 
atmospheric circulation) . 

4.3.2. The Effect of Dissipation on the RCB 

Ignoring dissipation, we obtained the limiting K XZfCt it 
as P c -> oo. With dissipation, K ZZtCI n occurs at finite P c , 
provided there is dissipation at the RCB. As noted above, 
the RCB occurs at 11 kbar with K zz — K ZZtCT n in Fig.[SJ 
Appendix iBl derives the relation between dissipation and 
the maximum P c in equation (IF33|) . 

We consider it physically desirable to restrict P c to 
finite pressures. Infinite P c obviously violates some of 
our idealizations, notably the plane parallel and ideal gas 
approximations. It is encouraging that dissipation alters 
P c without qualitatively changing the insights (namely 
K zz , CT it) gleaned from the non-dissipative model. 

Restricting P c to finite values was not crucial for the 
energetic balance arguments above. Though F c is larger 
for smaller P c , it is still too small to be the main factor 
that limits the eddy flux. 

To illustrate some of these points consider the case 
e = fbuoy, i-e. with no floor, e G , to the dissipation. In 
this case there is no dissipation at the RCB, and we find 
that eddy diffusion can still push P c — » oo. Nevertheless 
fbuoy by itself does still lower K ZZtCT it, by ~ 1/3. This 
can be seen in the e c — » limits (i.e. the left of the plot) 
of Fig. M 

4.4. Varying the Opacity Law and EOS 

All of our plots and numerical estimates have used an 
opacity law k oc P and an ideal gas EOS with V a d = 2/7. 
We provide general scalings for many results to show the 
effect of varying these parameters. 

Our qualitative results hold for all "reasonable" choices 
of power law opacities and EOS. As discussed in £12.11 
reasonable means that Voo > V a d and a > — 1 so that 
radiative equilibrium solutions become convectively un- 
stable at depth. We tested our results with an alternate 
opacity law n cx T 2 (a = and j3 — 2), that approxi- 
mates H~ opacities for T > 2000 K or dust grain opac- 
ities at colder temperatures. Like n oc P, this law also 
has = (1 + a)/(4 - 0) = 1/2 > V ad; 

The behavior of -ftT Z z, crit — shown in equation (|32[) 
— is particularly important for interpreting our results. 



The scaling with internal entropy is quite steep with 
K ZZyCr it oc Ti°- 5 and K zz ^- Xt oc Tj for our standard and 
alternate opacities, respectively. Generally, the entropy 
dependance becomes steeper for larger a (as in the above 
example) and also for a smaller V a d, but does not depend 
on p. Recall that the burial of the turbulent heat flux 
into the convective interior can bring an atmosphere with 
K zz > K zz>CI it towards energetic equilibrium. A steeper 
dependance of -K" Z z,crit on entropy means that less infla- 
tion is needed to enforce this equilbrium. 

The scaling of K gza a with Td CC p controls how our in- 
flation mechanism depends on the level of irradiation. 
Also the development of thermal inversions will lower 
Tdccp for a fixed level of irradiation. We find K ZZt „it oc 

T delp or K **,ait oc Tdccp' a § am for tne standard and 
alternate opacities, respectively. The Td ee p dependance 
becomes more steeply negative for larger a (again the 
dominant effect in our example) and also for larger /3 
(less important in our example) and smaller V a d- 

These scalings emphasize that a more detailed treat- 
ment of turbulent heat fluxes in hot Jupiters should in- 
clude realistic opacities and equations of state. Non- 
powerlaw be havior could have si gnificant consequences. 
For instance, IGuillot et al.l ([1994D demonstrated the im- 
portance of an opacity window near ~ 2000 K. This win- 
dow can give rise to an isolated convective layer sand- 
wiched between two radiative zones. It would be inter- 
esting to consider how a downward turbulent heat flux 
would interact with such a region. Future work that gen- 
eralizes our self-similar approach can address these more 
detailed issues. 

5. COMPARISON WITH PREVIOUS WORK 

5.1. Simulations of Atmospheric Circulation 

Hydrodynamic simulations of hot Jupiter atmo- 
spheres have been stu died in local ([Burkert et al.l [20051 : 
| Li fc Goodmanl 120101 herea f ter LG10) and glob a l (e.g . 
Cooper fc Showman! [20051: iRauscher fc Menqul [2010: 
IShowman et al.1120091: IDobbs-Dixon fc Lirj|2008f) models. 
A major goal of these studies is to determine the circu- 
lation induced by stellar irradiation. We first discuss es- 
timates of K zz from these simulations, and then address 
the question of whether radiatively forced turbulence ex- 
tends throughout the radiative zone, as we have assumed. 
Other sources of turbulence — perhaps involving mag- 
netic fields — may also exist, but we do not analyze them 
here. 

The global circulation models of IShowman et al.l 
(2009) estimate K zz ~ 10 11 cm 2 /s at mbar pressures. 
This does not contradict our constraint on upper atmo- 
spheric mixing from the efficiency of radiative forcing 
[equation (|33[) ] when extrapolated to such low pressures. 
Moreover at low pressures radiative losses can lower the 
eddy flux (see q3.1[) and further weaken our constraint 
on K zz . We caution that K zz estimates from GCMs are 
rough, since t hey only resolve relati vely large scale flow 
patterns. The lShowman et al.l ([20 09; K zz estimate arises 
from multiplying measured vertical speeds by the scale- 
height H . While a useful guide to what is possible, this 
does not constitute a direct measurement of diffusion. 

Local hydrodynamic simulations can better resolve 
turbulent flows, though as usual with Reynolds num- 
bers far lower than reality. The calculations of LG10 
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find an effective turbulent viscosity of v t ~ 0.001 - 

0. 01c s iJ - 10 10 — 10 11 cm 2 /s. It is unclear if this 
viscosity should be interpreted as a mixing coefficient. 
Their simulation with v t = 0.015c s -ff has an RMS ver- 
tical speed w ~ 0.3c s . Thus assuming v t ~ i^ zz is not 
consistent with the simple estimate K zz ~ u>£ because 
it gives a length scale ^ ~ H/20, smaller than the grid 
spacing of H/10. It should not be surprising that simple 
estimates based on 3D isotropic turbulence fail, given the 
organized structure in their 2D "turbulent" state. 

Moreover the forcing in LG10 was not by irradiation, 
but chosen to be large enough to drive super-sonic flows 
despite artificial viscosities that are large for numerical 
reasons. Thus attempting to interpret the Carnot effi- 
ciency of stellar irradiation may overextend their results. 
Despite these caveats, if the LG10 simulations apply near 
the mbar weather layer, there is again no contradiction 
with equation (|33j) . 

Global simulations indicate that the shear layer may 
not extend throughou t the radiative zone. For instance 
iShowman et al.l (|2009f ) find that strong (~km/s) zonal 
winds terminate at ~ 10 bar. They argue that circula- 
tion stops because the planet is horizontally isothermal 
at these depths, removing the local forcing. However, hot 
Jupiter atmospheres have a large separation between the 
radiative timescales in the weather layer and the deep 
radiative layer. It is possible (and arguably likely) that 
unresolved or long timescale dynamics could push the 
shear layer even deeper. This possibility should be ex- 
plored in future modeling studies. 

How far turbulence can extend below the shear layer 
is also uncertain. LG10 force a shear layer of with ~ 
2H, but find that turbulence (or at least some disordered 
motion) extends throughout their box of size ~ 5H (see 
their Fig. 10). Turbulence that extends a full 5H below 
a shear layer could thus extend to P ~ e 5 • 10 bar ~ 1.5 
kbar. These issues deserve further study, but it cannot 
be ruled out that the radiative zone is turbulent down to 
the RCB. 

5.2. Thermal Inversions via TiO Diffusion 

We now provide a more detailed interpretation of our 
results in terms of the S09 constraints on the mixing re- 
quired to loft TiO and create thermal inversions. We 
emphasize that the constraints on K ZZl and especially 
on the depth dependance of K zz depends sensitively on 
whether or not there is a cold trap. Wherever TiO con- 
denses, one must consider the mixing of dust grains, not 
just molecules. 

First consider the case where there is no cold trap, 
and TiO is always in the vapor phase. In the S09 analy- 
sis, this is only possible for the most intensely irradiated 
planet, WASP 12b, in part because thermal inversions 
lower Tdeep and favor condensation at depth. The mixing 
of TiO vapor requires K zz ~ 10 7 cm 2 /s at P ~ 1 mbar, 

1. e. the height of the inversion where TiO is needed. The 
specification of pressure level is important because the 
constraint on molecular diffusion scales as K zz oc P^ 1 in 
an isothermal atmosphere. Thus the constraint is weaker 
at depth. In a model where the actual K zz oc P^ 1 (as 
in t j4.2[) we only require K zz ~ 10 cm 2 /s at the kbar 
RCB. In this case our model predicts a minimal effect 
of the eddy flux on the planet's evolution, though the 



dissipation of turbulence above the RCB could still be 
significant. 

We briefly summarize how this K zz ~ 10 7 cm 2 /s limit 
and the depth dependance arises, and refer the reader 
to S09 for details. With only molecular viscosity, i.e. 
collisions, the TiO scaleheight would be hydrostatic, and 
thus ~ 30 times smaller than the dominant H2 species. 
S09 conclude that the turbulent diffusion must exceed 
the molecular diffusion Duo by a factor ~ 100, due to 
the many (~ 14) scaleheights between P ~ mbar and 
the kbar RCB. Estimating DxiO as usual — the product 
of the mean free path and thermal speed — gives an 
inverse scaling with density, and thus also with pressure 
in isothermal regions. The numerical value of Z?TiO ~ 
10 5 cm 2 /s at P ~ mbar, combined with the factor of 
100 excess needed for efficient turbulent mixing, gives 
the K zz ~ 10 7 cm 2 /s constraint. 

Now consider the case where grains do condense some- 
where in the radiative zone, which is true for most hot 
Jupitcrs (especially those with inversions). The K zz re- 
quired to loft TiO increases to 10 8 — 10 11 cm 2 /s depend- 
ing on the size of grains and the depth of the cold trap. 
A rough estimate of K zz ~ Vt eim L follows from equat- 
ing the diffusion timescale, L 2 /K zz , to the dust settling 
timescale L/vtcrun where L is the depth of the cold trap 
and Vterm is the grain's terminal speed. While the ter- 
minal speed will increase with particle size, it does not 
depend on atmospheric density for the relevant viscous 
(i.e. Stokes') drag. Thus unlike the case of molecular 
viscosity, the K zz required for mixing condensed grains 
does not decline with depth. 

We can thus conclude that the K zz required to mix 
TiO and create thermal inversions in a hot Jupiter are 
likely excessive, provided TiO condenses somewhere in 
the radiative zone. This follows from Fig. 2] which shows 
that K zz ~ 10 9 cm 2 /s appear off-scale. The internal 
(specific) entropy of a planet would have to increase by 
> 2/cb/^p for a planet with Td CC p ^ 1500 K, as is the 
case for all planets in S09. It seems likely that such large 
entropies would overinflate a planet, based on studies 
such at AB06. This is especially true because Fig. @] 
neglects dissipation which can further inflate a planet, 
which can occur by lowering K ZX>C[ & as shown in Fig. [8] 

However, we cannot firmly declare that the TiO hy- 
pothesis fails. This is largely due to the approximate 
nature of our treatment of opacities and the equation of 
state. A more conclusive analysis of the TiO hypothesis 
would require a detailed atmospheric model that includes 
the eddy fluxes and turbulent dissipation described in 
this work. Such a study would also have to abandon the 
fixed flux bottom boundary condition used in most at- 
mospheric models (including S09). Instead the bottom 
boundary should be a fixed adiabat, chosen to match 
the observed radius. This is subject as usual to assump- 
tions about composition and presence of a core. However 
allowing the flux to float is needed for a consistent deter- 
mination of the RCB. This is crucial for understanding 
the coupled relationship between compositional mixing 
and cooling history. 

6. CONCLUSIONS 
6.1. Summary of Results 
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We have investigated how forced turbulent mixing af- 
fects the energetic balance and structure of the stably 
stratified, radiative layers of hot Jupiters. It is crucial 
to understand how the radiative layer matches onto the 
convective interior at the radiative-convective boundary 
(RCB). This regulates the rate at which the planet cools 
and therefore controls the evolution of the planet's ra- 
dius. 

Previous work has invoked turbulent eddy diffusion of 
molecular species (and dust). This mixing changes the 
opacity to provide better- fitting model spectra of transit- 
ing planets, especially those that appear to have thermal 
inversions (S09). We find that turbulent mixing of this 
kind does not just redistribute chemical species but also 
significantly affects energetics. 

Forced turbulent mixing in stable, radiative regions 
drives a downward flux of energy that pushes the RCB 
deeper in the atmosphere, lowering the planet's cooling 
rate. We found an upper limit to the strength of turbu- 
lent diffusion, K Z z,crit> that can be achieved in steady- 
state. Beyond this limit, the downward flux of energy 
will heat the convective interior and inflate the planet. 
We did not directly model this entropy growth because 
our model was steady state and did not include over- 
shoot across the RCB. The deep, marginally stable layer 
in our solutions would not strongly inhibit overshoot, so 
heat burial by this mechanism is a likely outcome. Our 
solutions indicate that interior heating brings the planet 
back towards steady state, because higher entropy plan- 
ets have larger K ZZtCr i t (see Fig. |4j. Our mechanism is 
thus a "mechanical greenhouse effect" with the role of 
sunlight in the traditional greenhouse being played by 
forced turbulent mixing. 

For non-uniform turbulent mixing we find that our con- 
straint on K ZZ CI a applies near the RCB. More specifi- 
cally this constraint applies at a pressure Pdeep, which 
lies below the deep isothermal region where the radiative 
layer is transitioning towards convective instability. Our 
constraints on turbulence in the upper atmosphere are 
less stringent. However, the downward flux of mechani- 
cal energy likely cannot exceed a small fraction (probably 
at the percent level) of the stellar irradiation if it is sup- 
plied by weather-layer winds. If turbulence is too weak 
at the bottom of the radiative layer, it will not dredge up 
heavy molecular species — particularly those that con- 
dense onto dust grains — to serve as opacity sources near 
the observable photosphere. 

Turbulence also deposits heat in the radiative atmo- 
sphere when it decays. Non-turbulent sources of energy 
dissipation — including non-linear wave breaking and 
ohmic dissipation — also affect energetic balance. We 
find that including energy dissipation reduces K ZZiCI n, 
and thus makes it easier to inflate a planet for a given 
level of forced turbulence. 

We find a characteristic scale of K zz ^ CI i t ~ 10 3 - 
10 5 cm 2 /s for typical hot Jupiter parameters, even ignor- 
ing dissipation. This is orders of magnitude below values 
quoted in the literature of 10 7 — 10 11 cm 2 /s for the mix- 
ing o f chemical species (Spiegel ct al. 2009; iZahnle et al.l 
2009). We caution that our quantitative results should 
be taken as illustrative, due to the approximations de- 
tailed in jgl 

Thus when turbulence is important for the redistribu- 



tion of opacity sources, it also has significant energetic 
consequences. Our model represents a step towards more 
self-consistent modeling of exoplanet atmospheres. 

6.2. Applications and Extensions 

iGuillot fc Showma n (2002) first noted that dissipation 
in radiative layers can slow planetary contraction by 
pushing the RCB to higher pressures. They imagined 
dissipation concentrated in the upper regions of the at- 
mosphere, sourced by insolation driv en winds. Other au- 
thors (e.g. iBodenheimer et al.ll2003[ ) have included the 
IGuillot fc Showman (2002) dissipation prescription to 
model exoplanet radii. 

However we are unaware of previous works that in- 
clude the mechanical flux of energy due to turbulence or 
waves — our -F e ddy — in radiative layers. Incorporating 
this effect in detailed planetary evolution models would 
require a more precise treatment than this exploratory 
study. Notably it would require a global model with a 
realistic equation of state and opacity, not the power laws 
considered here. 

In this paper we treat diffusion and dissipation by 
forced turbulence as free parameters to allow a general 
analysis. This approach is justified since we currently 
lack a detailed understanding of these processes. Thus 
our model can be used to test how any specific model for 
turbulence and/or energy dissipation affects the struc- 
tu re and evolution of the plane t. For instance, the work 
of iBatygin fc Stevenson! (|2010h consider ohmic dissipa- 
tion in radiative and convective regions, but do not re- 
compute the effect of the dissipation on the structure of 
the radiative layer. Our point is not to critique, but to 
emphasize that a more consistent treatment of energetics 
may make it easier to inflate a planet — by any number 
of mechanisms. 

We also hope to include the effects of our study in de- 
tailed 3D radiation-hydrodynamical simulations of exo- 
planet atmospheres. This would involve adding sub-grid 
physical prescriptions for eddy fl uxes to global circulatio n 
models (GCM) such as th ose bv lShowman et all (|2009f k 
iRauscher fc Merioul (l2010fl and/o r non-hydrostatic s imu- 
lations (jDobbs-Dixon et afll2ojoh . iGoodmanl (|2009f ) dis- 
cusses the need to incl ude explicit sources of dissipation. 
iLi fc Goodman! (|2010D present a first step towards sub- 
grid modeling of turbulence due to Kelvin-Helmholtz in- 
stabilities. Breaking of vertically propagating gravity 
waves represent anoth er possible source of turbulence 
(Sh owman et al.ll2009| ). When all these effects are taken 
into account, the inflated radii of transiting planets may 
not be surprising after all. 

We have focused on hot Jupiters, but the physics we 
describe in principle applies to other atmospheres. Our 
limits on eddy diffusion become less severe as objects are 
more weakly irradiated (see Fig. |4|). Thus our mecha- 
nism would not inflate longer period gas giants, including 
Jupiter and Saturn. More intrinsically luminous objects 
like brown dwarfs will similarly be less affected. The 
expanding inventory of transiting exoplanets, at wider 
separations, is an excellent test of radius evolution mod- 
els. 

Finally we comment on a possible connection to 
the atmosphere of Venus. Venus has a marginally 
stable pseudo-adiaba t beneath a thick cloud deck 
(|Schubert et al.l [T980f ) . The clouds shield sunlight from 
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warming the surface of Venus. Because of this, Venus' at- 
mosphere would be nearly isothermal were it not for some 
poorly understood mechanical stirring process. Our well- 
stirred hot Jupiter atmospheres also have marginally sta- 
ble pseudo-adiabats at depth. Perhaps Venus is a very- 
well-stirred analog of a hot Jupiter. Obvious differences 
exist, for instance the non-negligible fraction of sunlight 
that reaches the Venusian surface versus the radiant flux 
escaping the core of hot Jupiters. Pursuing analogies 
such as these should improve our understanding of worlds 
near and far. 
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APPENDIX 

A. ZERO FLUX SOLUTIONS 

In addition to the solution technique described in §3.2[ we also used a more specialized technique that applies to a 
reduced model with no dissipation and a net flux F — 0. As we show in §4.11 this is an interesting case because the 
downward flux of heat can drive the core flux, F c — > 0, and the flux remains zero in the absence of dissipation. This 
alternate technique allows us to check our results. It allows integration from the top down, compared to the bottom 
up integration from the RCB (which recedes to infinite pressure for the zero flux solution). We do not need to specify 
the parameter ip c (eq. [26]), which diverges for these zero flux solutions. With the reduced model we can analytically 
explore solution properties. We will do this below to explain why we cannot obtain solutions with £ < —1.5 in §4.21 

Ignoring dissipation and setting F — we can rearrange equation (|23l) as 



dT T t V ad k lad T \- L 

We can integrate from P = with the boundary condition T(0) = Td CO p- The solutions approach an adiabatic profile 
T oc P Vad at large pressure (shown below). One could find solutions by integrating with various K zz values, iterating 
until you land on a desired adiabat. Instead we again take a self-similar approach, non-dimensionalizing the pressure 
using K zz and the opacity law. We choose a physical scale for the pressure by matching the self-similar solution onto 
a chosen adiabat. This fixes the value of K zz (at a reference pressure if K zz is not uniform). The special value of K zz 
that gives a zero flux solution for a given adia bat is the limiting K zz ^a discussed in £ 14.11 

We now show that the solution to equation (|A1[) becomes adiabatic at large pressure, at least if K zz is constant. If 
the final term in parenthesis vanishes we have V — > V a d- We can show that an adiabatic solution is consistent by 
assuming T oc pv» d and then confirming that the final term 

Y^L K ZJ! K p^Xvwv-D 1 (A2) 
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as P — > oo. With V'^ = (2+ a)/(5 — 0) we can show that the exponent in the last proportionality of equation (|A2|) is 
indeed negative. This is because VJ^ > Voo (which is generally true, for our opacity V x — 3/5 > 1/2) and because 
reasonable opacities have a > — 1 and Voo > V a d (see $2]). Thus a deep adiabat is a consistent solution at large 
pressure, for constant K zz . 
Now consider a spatially varying K zz oc P^. In this case the assumption of an adiabat at depth gives 



V 1 



V ad 1 + c c P< 



(A3) 



as P — ^ oo, where is constant and = V a d(5 — (3) — (2 + a + £)• Consistency (i.e. V — > V a d) requires that 6q < 
or ( > —(2 + a) + V a( j(5 — (3) = —11/7 w —1.57. More negative values of ( do not approach an adiabat even at infinite 
pressure. 

This explains the behavior of ( < —1.4 solutions in H4.2I However we should not over- interpret the physical 
significance of this result since it may merely reflect a limitation of the self-similar approach. More study is needed 
of what happens when a large mechanical flux of energy enters the top of the atmosphere, but turbulent heat burial 
is very weak near the RCB (which is the case for negative (). For the time being, we thus consider the irradiation 
efficiency constraint in equation (|33j) as our best limit on turbulent mixing in the upper atmosphere. 

B. MAXIMUM DISSIPATION AT RCB 

This appendix derives how dissipation at the RCB restricts the RCB depth P c . We use the consistency requirement 
that dV/dP > at the RCB. This states that the solution will indeed be convectively stable above the RCB. We 
rearrange equation (|2"31 to give 

_V_ = F + F iso 

V ad V ad fc rad T/P + F iso 1 ; 

and at the RCB (denoted by the 'c' subscript) we have F = F c = V a d(fcradP/P)c- The gradient of V at the RCB is 



dP \yZ) c ~ F c + P iso , c dP 



V a dfc r adP 



P 



(B2) 



where the derivatives of Fi so cancel. Thus the limit on dissipation at the RCB, using equation (|18l) and requiring 
(dV/dP) c > 0, is 

ec - (1 + a) l 1_ v^J^ = W (B3) 
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bmce F c /P c oc P c ' , the maximum depth of the RCB declines with increasing dissipation at the RCB. 

This result is very useful in finding -ftT 2ZiC rit values with dissipation. Without dissipation we could easily find A^crit 
by increasing ip c (eq. [26]) to arbitrarily large values (which results in Fig. [3]), or by using the zero flux solutions of 
appendix [3] Instead we a dd a twist to the self similar technique of £13.21 We set the dissipation at the RCB to the 
limiting value of equation (|B3|) . We further assume the dissipation is of the form of equation (f30]) so that the strength 
of the dissipation is set by the dimensionless / e . We have thus coupled dissipation to K zz mathematically. (It doesn't 
matter if they are unrelated physical ly, s ince f e can be adjusted to give any desired level of dissipation.) Then we can 
solve equations (gU), O, ® and (|B3| for 

V 7/ £ V ad T c V ' 

Since T c is not known until we obtain a solution, we use an iteration procedure: guess a value for T c (but start with 
something too low), use the estimate of ip c from equation (IB4[) to integrate the model equations (f2"3")l and (|24p . use the 
resulting T c to refine t/'c and repeat. Though a bit convoluted, this procedure converges. 

Finally note that our co nstr aint on dissipation at the RCB does not ensure stability at all P < P c . Assume that 
the inequality in equation (IB3|) holds so that there is a radiative layer above the RCB. Dissipation could still create 
another convective layer at greater height. However our self-similar solutions cannot include such structures. Thus we 
leave it to a future work to consider these more complicated scenarios and their effects on planetary evolution. 



